The use of hybrid cellular automaton models for improving cancer therapy, In Proceedings, Cellular Automata: 6th International Conference on Cellular Automata for Research and Industry, ACRI 2004, Amsterdam, The Netherlands, eds P.M.A. Sloot, B. Chopard, A.G. Hoekstra

The Hybrid Cellular Automata (HCA) modelling framework can be an efficient approach to a number of biological problems, particularly those which involve the integration of multiple spatial and temporal scales. As such, HCA may become a key modelling tool in the development of the so-called intergrative biology. In this paper, we first discuss HCA on a general level and then present results obtained when this approach was implemented in cancer research

A mathematical model of Doxorubicin treatment efficacy on non-Hodgkin’s lymphoma: Investigation of current protocol through theoretical modelling results

Doxorubicin treatment outcomes for non-Hodgkin’s lymphomas (NHL) are mathematically modelled and computationally analyzed. The NHL model includes a tumor structure incorporating mature and immature vessels, vascular structural adaptation and NHL cell-cycle kinetics in addition to Doxorubicin pharmacokinetics (PK) and pharmacodynamics (PD). Simulations provide qualitative estimations of the effect of Doxorubicin on high-grade (HG), intermediate-grade (IG) and low-grade (LG) NHL. Simulation results imply that if the interval between successive drug applications is prolonged beyond a certain point, treatment will be inefficient due to effects caused by heterogeneous blood flow in the system

An Integrated Disease/Pharmacokinetic/Pharmacodynamic Model Suggests Improved Interleukin-21 Regimens Validated Prospectively for Mouse Solid Cancers

Interleukin (IL)-21 is an attractive antitumor agent with potent immunomodulatory functions. Yet thus far, the cytokine has yielded only partial responses in solid cancer patients, and conditions for beneficial IL-21 immunotherapy remain elusive. The current work aims to identify clinically-relevant IL-21 regimens with enhanced efficacy, based on mathematical modeling of long-term antitumor responses. For this purpose, pharmacokinetic (PK) and pharmacodynamic (PD) data were acquired from a preclinical study applying systemic IL-21 therapy in murine solid cancers. We developed an integrated disease/PK/PD model for the IL-21 anticancer response, and calibrated it using selected “training” data. The accuracy of the model was verified retrospectively under diverse IL-21 treatment settings, by comparing its predictions to independent “validation” data in melanoma and renal cell carcinoma-challenged mice (R2>0.90). Simulations of the verified model surfaced important therapeutic insights: (1) Fractionating the standard daily regimen (50 µg/dose) into a twice daily schedule (25 µg/dose) is advantageous, yielding a significantly lower tumor mass (45% decrease); (2) A low-dose (12 µg/day) regimen exerts a response similar to that obtained under the 50 µg/day treatment, suggestive of an equally efficacious dose with potentially reduced toxicity. Subsequent experiments in melanoma-bearing mice corroborated both of these predictions with high precision (R2>0.89), thus validating the model also prospectively in vivo. Thus, the confirmed PK/PD model rationalizes IL-21 therapy, and pinpoints improved clinically-feasible treatment schedules. Our analysis demonstrates the value of employing mathematical modeling and in silico-guided design of solid tumor immunotherapy in the clinic

Predicting Outcomes of Prostate Cancer Immunotherapy by Personalized Mathematical Models

Therapeutic vaccination against disseminated prostate cancer (PCa) is partially effective in some PCa patients. We hypothesized that the efficacy of treatment will be enhanced by individualized vaccination regimens tailored by simple mathematical models.We developed a general mathematical model encompassing the basic interactions of a vaccine, immune system and PCa cells, and validated it by the results of a clinical trial testing an allogeneic PCa whole-cell vaccine. For model validation in the absence of any other pertinent marker, we used the clinically measured changes in prostate-specific antigen (PSA) levels as a correlate of tumor burden. Up to 26 PSA levels measured per patient were divided into each patient's training set and his validation set. The training set, used for model personalization, contained the patient's initial sequence of PSA levels; the validation set contained his subsequent PSA data points. Personalized models were simulated to predict changes in tumor burden and PSA levels and predictions were compared to the validation set. The model accurately predicted PSA levels over the entire measured period in 12 of the 15 vaccination-responsive patients (the coefficient of determination between the predicted and observed PSA values was R(2) = 0.972). The model could not account for the inconsistent changes in PSA levels in 3 of the 15 responsive patients at the end of treatment. Each validated personalized model was simulated under many hypothetical immunotherapy protocols to suggest alternative vaccination regimens. Personalized regimens predicted to enhance the effects of therapy differed among the patients.Using a few initial measurements, we constructed robust patient-specific models of PCa immunotherapy, which were retrospectively validated by clinical trial results. Our results emphasize the potential value and feasibility of individualized model-suggested immunotherapy protocols

The Local Optima Level in Chemotherapy Schedule Optimisation

In this paper a multi-drug Chemotherapy Schedule Optimisation Problem (CSOP) is subject to Local Optima Network (LON) analysis. LONs capture global patterns in fitness landscapes. CSOPs have not previously been subject to fitness landscape analysis. We fill this gap: LONs are constructed and studied for meaningful structure. The CSOP formulation presents novel challenges and questions for the LON model because there are infeasible regions in the fitness landscape and an unknown global optimum; it also brings a topic from healthcare to LON analysis. Two LON Construction algorithms are proposed for sampling CSOP fitness landscapes: a Markov-Chain Construction Algorithm and a Hybrid Construction Algorithm. The results provide new insight into LONs of highly-constrained spaces, and into the proficiency of search operators on the CSOP. Iterated Local Search and Memetic Search, which are the foundations for the LON algorithms, are found to markedly out-perform a Genetic Algorithm from the literature

Dickkopf1 Regulates Fate Decision and Drives Breast Cancer Stem Cells to Differentiation: An Experimentally Supported Mathematical Model

BACKGROUND: Modulation of cellular signaling pathways can change the replication/differentiation balance in cancer stem cells (CSCs), thus affecting tumor growth and recurrence. Analysis of a simple, experimentally verified, mathematical model suggests that this balance is maintained by quorum sensing (QS). METHODOLOGY/PRINCIPAL FINDINGS: To explore the mechanism by which putative QS cellular signals in mammary stem cells (SCs) may regulate SC fate decisions, we developed a multi-scale mathematical model, integrating extra-cellular and intra-cellular signal transduction within the mammary tissue dynamics. Preliminary model analysis of the single cell dynamics indicated that Dickkopf1 (Dkk1), a protein known to negatively regulate the Wnt pathway, can serve as anti-proliferation and pro-maturation signal to the cell. Simulations of the multi-scale tissue model suggested that Dkk1 may be a QS factor, regulating SC density on the level of the whole tissue: relatively low levels of exogenously applied Dkk1 have little effect on SC numbers, whereas high levels drive SCs into differentiation. To verify these model predictions, we treated the MCF-7 cell line and primary breast cancer (BC) cells from 3 patient samples with different concentrations and dosing regimens of Dkk1, and evaluated subsequent formation of mammospheres (MS) and the mammary SC marker CD44(+)CD24(lo). As predicted by the model, low concentrations of Dkk1 had no effect on primary BC cells, or even increased MS formation among MCF-7 cells, whereas high Dkk1 concentrations decreased MS formation among both primary BC cells and MCF-7 cells. CONCLUSIONS/SIGNIFICANCE: Our study suggests that Dkk1 treatment may be more robust than other methods for eliminating CSCs, as it challenges a general cellular homeostasis mechanism, namely, fate decision by QS. The study also suggests that low dose Dkk1 administration may be counterproductive; we showed experimentally that in some cases it can stimulate CSC proliferation, although this needs validating in vivo

From regional pulse vaccination to global disease eradication: insights from a mathematical model of Poliomyelitis

Mass-vaccination campaigns are an important strategy in the global fight against poliomyelitis and measles. The large-scale logistics required for these mass immunisation campaigns magnifies the need for research into the effectiveness and optimal deployment of pulse vaccination. In order to better understand this control strategy, we propose a mathematical model accounting for the disease dynamics in connected regions, incorporating seasonality, environmental reservoirs and independent periodic pulse vaccination schedules in each region. The effective reproduction number, $R_e$, is defined and proved to be a global threshold for persistence of the disease. Analytical and numerical calculations show the importance of synchronising the pulse vaccinations in connected regions and the timing of the pulses with respect to the pathogen circulation seasonality. Our results indicate that it may be crucial for mass-vaccination programs, such as national immunisation days, to be synchronised across different regions. In addition, simulations show that a migration imbalance can increase $R_e$ and alter how pulse vaccination should be optimally distributed among the patches, similar to results found with constant-rate vaccination. Furthermore, contrary to the case of constant-rate vaccination, the fraction of environmental transmission affects the value of $R_e$ when pulse vaccination is present.Comment: Added section 6.1, made other revisions, changed titl

Patterns in Age-Seroprevalence Consistent with Acquired Immunity against Trypanosoma brucei in Serengeti Lions

Trypanosomes cause disease in humans and livestock throughout sub-Saharan Africa. Although various species show evidence of clinical tolerance to trypanosomes, until now there has been no evidence of acquired immunity to natural infections. We discovered a distinct peak and decrease in age prevalence of T. brucei s.l. infection in wild African lions that is consistent with being driven by an exposure-dependent increase in cross-immunity following infections with the more genetically diverse species, T. congolense sensu latu. The causative agent of human sleeping sickness, T. brucei rhodesiense, disappears by 6 years of age apparently in response to cross-immunity from other trypanosomes, including the non-pathogenic subspecies, T. brucei brucei. These findings may suggest novel pathways for vaccinations against trypanosomiasis despite the notoriously complex antigenic surface proteins in these parasites